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ABSTRACT 



Context. Stellar activity data provide evidence of activity wave branches propagating polewards rather than equator- 
wards (the solar case). This evidence is especially pronounced for the well-observed subgiant HR 1099. Stellar dynamo 
theory allows polewards propagating dynamo waves for certain governing parameters. 
Aims. We try to unite observations and theory. 

Methods. Taking into account the preliminary stage both of observations of polar activity branches and of the determi- 
nation of the governing parameters for stellar dynamos, we restrict our investigation to the simplest mean-field dynamo 
models, while recognizing more modern approaches to be an essential development. 

Results. We suggest a crude preliminary systematization of the reported cases of polar activity branches. Then we 
present results of dynamo model simulations which contain magnetic structures with polar dynamo waves, and identify 
the models which look most promising for explaining the latitudinal distribution of spots in dwarf stars. Those models 
require specific features of stellar rotation laws, and so observations of polar activity branches may constrain internal 
stellar rotation. Specifically, we find it unlikely that a pronounced poleward branch can be associated with a solar-like 
internal rotation profile, while it can be more readily reproduced in the case of a cylindrical rotation law appropriate 
for fast rotators. We stress the case of the subgiant component of the active close binary HR 1099 which, being best 
investigated, presents the most severe problems for a dynamo interpretation. Our best model requires dynamo action 
in two layers separated in radius. This interpretation requires some change in the paradigm of stellar magnetic studies, 
as it explains surface manifestations in a subgiant as a joint effect of shallow and deep layers of the stellar convective 
zone, rather than of a surface magnetic field only, as appears to be the case in dwarf stars. 

Conclusions. Observations of polar activity branches provide valuable information for understanding stellar activity 
mechanisms and internal rotation, and thus deserve intensive observational and theoretical investigation. Current stel- 
lar dynamo theory seems sufficiently robust to accommodate the phenomenology. 

Key words. Stars: activity - Stars: magnetic field - Dynamo - Sun: activity - Sun: magnetic topology 



1. Introduction 

It is widely accepted that the solar activity cycle is more 
than just a quasiperiodic variation of sunspot number, be- 
ing rather an activity wave that propagates from mid so- 
lar latitudes towards the solar equator. A solar dynamo 
based on the joint effects of differential rotation and mirror- 
asymmetric convective motions in the form of the so-called 
a-effect (possibly together with meridional circulation) is 
considered to be the underlying mechanism for the prop- 
agation of activity waves. Indeed, this mechanism gives 
an equatorwards propagating wave of large-scale magnetic 
field for a suitable choice of the parameters governing dy- 
namo action. It is natural to expect that such a phe- 
nomenon will appear in a variety of stars with convective 
envelopes, and we might thus be led to expect equatorwards 
waves of stellar activity. 

In fact cycHc activity is k nown now for many 
stars of various sp ectral types, e.g. Baliunas et al. (1995)[ 
Olah et al. (2009) Clarification of the spatial configura- 
tion of the assumed activity wave is a much more deli- 



cate undertaking. However contemporary astronomy pos- 
sesses a range of tools, such as the technique of Doppler 
Imaging (hereafter DI), with which to investigate the 
problem. A comprehensive investigation of the problem 
still remains a desirable milestone for stellar astronomy; 
however some early results are already available , e.g. 



Berdyugina fc Henry (2007) Katsova et al. (2010) The 



Send offprint requests to: D.Moss 



point here is that at least some stars exhibit an activity 
wave that propagates polewards. 

For instance, the K-type subgiant component of the 
RS CVn system HR 1099 has been extensively studied 
through DI by, e.g., Vogt et al. (1999) and shows migration 
of spots from mid-latitudes towards tne rotational poles on 
a timescale of a few years. Indirect evidence of the same 
phenomenon is found for several late-type main-sequence 
stars and young solar analogues, from chromospheric line 
flux monitoring or photometric optical monitoring, respec- 
tively (see Sect. [2|). 

In mean- field dynamo models, the direction of migration 
of the large-scale magnetic field features depends in prin- 
ciple on two key factors - the sign of the a coefficient in 
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the relevant hemisphere and the radial gradient of angular 
velocity. 

The situation is however not so straightforward. The 
point is that in addition to the equatorwards branch demon- 
strated by sunspots, the solar activity displays a relatively 
weak poleward s branch, seen in some other tracers, e.g. 
polar faculae - Makarov & Sivaraman (1989) It looks im- 
plausible a priori that such weak additional branches could 
be responsible for stellar polewards branches, however the 
possibility should be recognized. 

The aim of this paper is to investigate how the re- 
cent and current observational data concerning polewards 
branches of stellar activity might be connected with ideas 
from stellar dynamo theory. We appreciate that the obser- 
vational situation after these pioneering results still remains 
quite uncertain, and so we study just the most traditional 
forms of stellar dynamos, i.e. mean- field dynamos based 
on differential rotation and a-effect with simple algebraic 
quenching. 

More recent ideas in solar dynamo theory, such as flux 
transport dyna mo s based on meridional circulation, e.g. 
Durney (1995)1 |Choudhuri, Schussler fc Dikpati (1995)[ 



Uikpati fc Gilman (2005J 
oT dynamo 



or 



saturation, e.g. 



dynamical schemes 
jKleeorin et al. (2003)[ 
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Subramanian fc Brandenburg (2004)"] are certainly 
likely to be important. We believe however that a simple 
initial approach is desirable and therefore we consider a 
classic dynamo wave model with a simple nonlinearity as 
the basic model in our research. In this model, the activity 
wave propagation is primarily associated with the joint 
action of differential rotation and a-effect rather than any 
effects of meridional circulation. It appears however that 
some observational features of polewards propagation are 
difficult to reproduce with this simple model. Therefore, 
we also investigate briefly some effects of meridional 
circulation. 

We also note the possible role of low order dynamo 
models in el ucidating some aspects of ste llar magnetic be- 
haviour, e.g. Wilmot-Smith et al. (2006) but here we con- 
centrate on models that we feel are more directly inter- 
pretable physically. 

Previous theoretical investigations of stellar dynamos 
have focussed on reproducing the dependence of activity cy- 
cle periods on ste llar parameters, in particular the rotation 
period (see, e.g., |Noyes et al. (1984)] Ossendrijver (1997)] 
Jouve et al. (2010)| and references therein) or on explaining 
high-latitu de or polar spots that are not observed in the 
Granzer et al. (2000) [ [Holzwarth et al. (2006)| 



Sun 



Igik et al. (2007)[ )). We present here for the first time a ten- 
tative systematization of the stellar butterfly diagrams that 
are now emerging from the observations, and try to explain 
the different behaviours by means of a simple mean-field 
dynamo model. 

We stress that our primary aim is not to produce 'defini- 
tive' mean-field models for any of the observed behaviours. 
Rather we attempt to illustrate the degree of uncertainty 
inherent in mean-field parametrizations, and also to show 
that many behaviours are reproducible by such models. 
Whilst the fundamental shortcomings of mean-field theory 
have attracted much interest, rather less attention has been 
given to, for example, investigating differences in behaviour 
caused by modest changes to parametrizations. This is, of 
course, two-edged. It reduces any predictive power of mean- 



field modelling, but also illustrates the possibility of ex- 
plaining non-standard behaviours by more unusual regimes. 

2. Butterfly diagrams and polewards activity waves 
for solar-like stars 

We introduce some tentative systematization of the infor- 
mation concerning the migration of activity patterns as de- 
rived from available observations. Distinct cases are sum- 
marized in Table (TJ which lists also our dynamo models 
with features resembling those observed (see Sect. 2]), and 
are briefly described below. They are listed in order of in- 
creasing Rossby number (i?o), that is the ratio of the rota- 
tion period of the star to the convective turnover time at 
the base of the convection zone; Ro can be roughly related 
to t he dynamo number in the bulk of the convection zone 
(cf. [Noyes et al. (1984)D . 



Case I: Among the RS CVn systems, the K-type compo- 
nent of HR 1099 has a record of DI maps extending over 
about twenty years with simultane ous coverage in wide - 
band optical photometry (cf. , e.g., Strassmeier (2009)). 
Berdyugina fc Henry (2007)] extending previous work 
by |Lanzaet al. (2006)[ built maps of the distribution 



of starspots on the active K-type subgiant. Two main 
active regions were found, one migrating from high 
latitudes (« 70°) towards mid-latitudes (w 40°), the 
other from mid-latitudes (« 40°) towards high latitudes 
(« 70°), these occurring more-or-less simultaneously. 
Several other RS CVn binaries show a general be- 
haviour similar to that of HR 1099, although their 
DI and photometric time series are less extended or 
have a more limited simultaneous coverage. A char- 
acteristic of the active components of RS CVn bi- 
nary systems is the presence of a polar spot that per- 
sists with little modification over timescales of decades, 
and which may be due to the advection of mag- 
netic flux to the polar regi on of the star by diffu- 



sion or meridional flows, e.g. Schrijver fc Title (2001) 
Mackay et al. (2004)[ [Holzwarth et al. (2006)[ In con 



trast, starspots at intermediate and low latitudes seem 
to follow a cyclic migration t hat might be associ ated 



with an oscillating dynamo, cf. Strassmeier (2009) 
Case II: In young, rapidly rotating stars, such as AB Dor 
(P,„t = 0.51 days) and LQ Hya (P^ot = 1-66 days), DI 
has revealed the simultaneous presence of spots at high 
and low latitudes, as well as a polar spot which has not 
been observed in all seasons, thus indicating a less per- 
siste nt feature than in the case of the RS CVn systems 
(e.g., Kovari et al. (2004) ). In AB Dor, the spots at low 
and inter mediate latitudes do not appear to migrate sig - 



niflcantly ( [Jarvinen et a l. (2005)[ [Jeffers et al. (2007)] , 
in contrast to the case of HR 1099 (case 1 above). 
Case HI: Solar-like stars with a rotation period of about 
5 — 40 days have been studied through the long- 
term monitoring of their chromospheric flux vari- 
ations, mainly in the framework of the classic 
Mt. Wilson H&K project and its recent exten- 



sions, e.g. Baliunas et al. (1995) [Baliunas et al. (1998 



Hall fc Lockwood (2004)[ |Donahue fc Baliunas 

report that 36 stars out of about 100 have several deter- 
minations of their rotation period extending over sev- 
eral seasons. Among them, 21 show patterns of rotation 
that vary with time or with the phase of the activity 
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cycle. Specifically, 12 stars display a pattern that re- 
sembles what would be expected from the solar butter- 
fly diagram, although in six of them the rotation period 
increases as the cycle progresses, in contrast to the so- 
lar case. One of the best examples is HD 114710, see 



Donahue & Baliunas (1992) 

'I'he dwarf stars showing anti-solar behaviour seem un- 
likely to possess an anti-solar pattern of surface differ- 
ential rotation, i.e. with the poles rotating faster than 
the equator, because this has never been found from DI 
observations ( Barnes et al. (2005)[), or from theoretica l 
models of stellar rotation (e.g. Riidiger et al. (1998)). 
Therefore, a plausible explanation is that their active 
regions migrate polewards rather than equator wards, 
which is what we deflne as case III. 
Case IV: There is another possibility to produce the phe- 
nomenology described in Case III. If stellar activity is 
not confined to latitudes close to the equator, but is 
well extended towards the poles, in addition to a pole- 
wards dynamo wave there can be another wave propa- 
gating from intermediate latitudes towards the equator. 
Which of the two waves dominates the modulation of 
the stellar flux depends on the inclination of the rota- 
tion axis with respect to the line of sight. If the star 
is viewed approximately pole-on, the polewards branch 
will dominate and the observed behaviour is anti-solar, 
while if the inclination is low, the star shows a solar-like 
behaviour because the low-latitude branch dominates, 
as suggested by, e.g. 



Ca ses V and VI: Among 



Messina k Guinan (2003) 
tEe 



stars considered by 
there are four stars 
trend of rotation period 



Donahue fc Baliunas (1994) 
that seem to reverse the 
variation at mid-cycle; six stars that show two narrow, 
but well separated bands of rotation, suggesting two 
stationary active latitude belts - we define this as case 
V; and, finally, stars that have hybrid patterns with 
one band showing a variation of the rotation period 
versus the cycle phase, while the other remains fairly 
constant - we take this behaviour to be representative 
of case VI. Stars with one or two fixed rotation periods 
as determined from Ca II H&K monitoring might 
be characterized by a standing dynamo wave, e.g. 



Baliunas et al. (2006) 



Cases III, IV, V, and VI could be different manifes- 
tations of the same kind of activity behaviour, which ap- 
pear distinct, either because of a different inclination of the 
stellar rotation axis which emphasizes either the polar or 
the equatorial region of a star, and/or a different inten- 
sity of the activity, or a phase shift between the polewards 
and the equatorwards branches of the butterfly diagram 
during the activity cycle. The available observations are 
still too limited to arrive at any sound conclusion on this 
point. Therefore, we prefer to consider all the suggested 
behaviours separately because each case makes a specific 
requirement for the theoretical butterfly diagram. 

3. Mean-field dynamos 

Now we address the problem from the other aspect 
and discuss how the butterfly diagrams appear in an 
assortment of stellar and solar-like dynamo models. 
We discuss here the simplest cases from the view- 
point of dynamo theory, i.e. standard mean-field dy- 



namos, e.g. Riidiger & Hollerbach (2004) with conven- 
tio nal boundary conditions and numerical implementation, 
e.g. Brandenburg & Subramanian (2005) 

We investigate solutions of the standard mean field dy- 
namo equation: 



dB 

'dt 



V X (u X B + aB - rjV X B), 



(1) 



where rj is the turbulent diffusivity and a represents the 
usual isotropic alpha-effect. The velocity field u = rtzu<j> + 
Mm, where is the angular velocity of rotation, zu the dis- 
tance from the rotation axis, (f) the unit vector in the az- 
imuthal direction, and itm the meridional flow. We restrict 
our investigation to axisymmetric solutions and solve the 
dynamo problem in a spherical shell, tq < r < 1, where r 
is the fractional radius. We adopt tq = 0.64 for many of 
the models, but we also investigate models with a deeper 
dynamo region with tq — 0.2. 

When modelling stellar magnetic fields, it is necessary 
to ensure that the field in the interior joins smoothly on to 
a force-free field in the external, very low density, region. 
Splitting the magnetic field into poloidal and toroidal parts, 
B = Bp + Bt, the Lorentz force can be written as 



L = (V X Bp) X Bp + (V X Bt) x Bt 
+ (V X Bt) X Bp + (V X Bp) x Bt. 



(2) 



For an axisymmetric field, the last term is identically zero, 
the first two are poloidal vectors and the third is toroidal. 
In the present case Bt = B^,. The condition L = can 
be satisfied by setting B^ — and V x Bp = 0; this pro- 
vides the boundary condition on the interior field applied at 
r = 1. At the lower boundary we use 'overshoot' boundary 
conditions, simulating the decay of the field to zero over a 
skin depth 6, in the form dg/dr = g/S, where g represents 
the azimuthal component of the vector potential for the 
poloidal field, or the toroidal field. These boundary condi- 
tions have been used before, and have been shown to have 
no significant effect on the results, except to reduce some 
field gradients near the base of the convective region. 

Models with a solar-like rotation law, based on that de- 
rived from helioseismology, and also models with a quasi- 
cylindrical rotation law appropriate to rapidly rotating 
lower mass dwarfs, are studied (Fig. [T]), with a variety of 
choices for a(r, 0) = Cafi{r)f2{9)/[l+{B/Bof], i.e. a naive 
a-quenching nonlinearity is used, where 9 is the colatitude 
measured from the North pole and Bq a reference magnetic 
field (see below). 

Different physical mechanisms can co-operate to pro- 
duce the a-effect and theoretical or observational (e.g. 
|Zhang et al. 2010] ) knowledge of its radial and colatitu- 
dinal distributions, here expressed through the functions 
/i(r) and /2(^), is quite preliminary. Therefore, we adopt 
only simple parametrizations and explore a number of op- 
tions to investigate the sensitivity of butterfly diagrams 
to the underlying assumptions. The turbulent diffusivity 
7/ is uniform in the outer part of the convection zone 
(r > 0.8), but decreases linearly to one half that value 
in the domain r < 0.7. Below the CZ proper we expect 
an overshoot region, where the turbulent intensity, corre- 
spondin g to the turb ulent resistivity, is further reduced. 
In Moss fc Brooke (2000) (which uses the Malkus-Proctor 
feedback onto the diti'erential rotation as the sole nonlin- 
earity, as opposed to the algebraic alpha-quenching here). 
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Table 1. A systematization of the cases where polewards migration or other features different from the solar case are 
deduced to explain observations. EW means equatorwards propagation, PW polewards propagation. If our investigations 
result in a model that appears suitable to explain a case, we identify it in the third column. 



Case 


Observational features/example 


Dynamo interpretation 


I 


PW migration from mid-latitude and EW migration 
from high latitude (HR 1099) 


See Sect. 13.2.51 


II 


Spot pattern extended in latitude 
with high and low latitude spots, 
but no definite migration during the cycle (AB Dor) 


Model 24 


III 


PW migration, i.e. the spot rotation period increases 
as the cycle progresses, contrary to the solar case 


Model 13c 


IV 


EW at low latitudes, PW at high latitudes 

(the latter can dominate when the star is viewed pole-on) 


Model 8 


V 


Two separated narrow activity bands 


Possible after a further 

specialization of stellar hydrodynamics 


VI 


Migration + a standing pattern 


Model 19 



only token recognition of this effect was made for computa- 
tional reasons. We have followed the same procedure here, 
recognizing that the gradient in turbulent diffusivity should 
be larger. Another way of regarding this is to consider the 
model as having a rather deeper CZ, extending to radius 
tq. From this viewpoint, the lower boundary condition, that 
allows a penetration of the field with a skin depth of the 
order of S, corresponds to the effect of a strongly reduced 
diffusivity immediately below the boundary. 

We make the dynamo equation non-dimensional in 
terms of the stellar radius R, the diffusion time 
R^/rjQ, where t^q is the maximum turbulent diffusiv- 
ity, and the mag netic field defined as in Sect. 3 of 
Moss & Brooke (2000) Thus we introduce the standard dy- 
namo numbers, C'q = aoR/r/o,Cui — floR^/rjo, where ao is 
a typical value of a and Qq is the maximum value of the an- 
gular velocity. In the ao; approximation, we can define the 
combined dynamo number D — CaC'i^, and this remains a 
useful quantity even in a^ui models. 

To integ rate the dynamo equat ion, we use the code 
described in Moss & Brooke (2000) which uses a Runge- 
Kutta integrator over a standard mesh with 61 points over 
ro < r < 1, and 101 points over < 9 < n, equally spaced. 
The results are described in Sect. 13.21 

A useful simplification of the general axisymmetric 
mean- field dynamo equation ([1]), known as the Parker mi- 
gratory dynamo, is considered in Sect. [XT] and is written in 
a standard non-dimensional form as: 



(3) 



(4) 



where / is a factor that allows for a latitudinal variation 
of the radially averaged angular velocity (see Sect. 13. ip . a 
is a nondimensional measure of the a effect, D = CaC^ 
is the dynamo number, B denotes the toroidal magnetic 
field and A the toroidal component of the vector poten- 
tial for the poloidal field. Both latter quantities are av- 
eraged in radial direc tion over the convective shell, see 



Baliunas et al. (2006) In other words, here the explicit ra- 



diai dependence has been removed, and the terms involving 
fi^ represent radial diffusion in a spherical shell of thick- 
ness approximately of the outer radius of the shell - 



e.g. /i w 3 is appropriate for the solar convection zone. We 
solve Eqs. ^ and (g]) in the domain < 9 < tt with A = 0, 
B = at the boundaries. 



3.1. Polar branches in the ID Parker dynamo 

We begin with a simple cartoon which explains the idea of 
stellar dynamos, i.e. with the ID Parker (1955) dynamo. 

First of all, we estimate the latitudinal variation of il, 
averaged over fractional radii (0.69,1) from a realistic solar 
rotation law (Fig. 1, left hand panel). Normalized to the po- 
lar value of the gradient, the modulation compared to that 
with no latitudinal variation is modelled by a function f{9) 
(Fig. [5]). Then we ran the Parker dynamo with the usual 
dynamo number, D = Dq say, replaced by _D = Dof{9), so 
f{9) = 1 gives the standard case. 

We take a = cos 6* sin™ 6*, with m = 0,2,4. The 
larger the value of m, the more concentrated around the 
equator the a effect. This simple para metrization has 
been 



used by other authors, e.g. |Riidiger et al. (2003) 

Charbonneau (2010)[ and we also adopt it here. We stress 
the inherent uncertainty in the spatial depend ence of a and 



refer to, e.g., Riidiger & Brandenburg (1995) for some jus- 



tification in the framework of mean held theory with a 
specific turbulence model. Specifically, they explore an a 
effect dependence with m = 2 that is suggested by the ex- 
tension of their turbulence theory to third order terms in 
fl • C/ oc cos 9, where is the angular velocity vector of 
the star and U the vector of the gradient of the turbulent 
diffusivity that points in the radial direction. In principle a 
further extension of the theory to the fifth degree in ft ■ U 
would introduce terms proportional to sin** 9 as we assume 
in our simple parametrization with to = 4 (of course we 
should then also consider expressions for alpha containing 
a combination of these dependences, but we decided that 
this was beyond the scope of this paper). Note also that 
the a effect may have a component arising from magne- 
tostrophic waves excited in the layers where the toroidal 
field is amplified and stored. This would lead to terms with 
TO > in the parametriz ation of the latit udinal dependence 
of the a effect (see, e.g., Schmitt (1987)). 



We emphasize that calculations of a from turbulence 
models are necessarily severely truncated, and stellar dy- 
namos operate in regimes remote from those in which such 
calculations are valid. Also, if we consider the potentially 
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Fig. 1. Equally spaced isorotation contours for the solar 
rotation law used (left hand panel) and contours for our 
quasi-cylindrical law for rapid rotators (right hand panel). 




Fig. 2. Function /(=< ^ >) used in the Parker model 
(Sect. EH). 



more useful, but more problematical, method of obtaining 
parametrizations of a from analysis of direct numerical sim- 
ulations, it would be surprising if m = or m = 2 or a com- 
bination was adequate. However such determinations are 
currently contentious, in spite of substantial progress, from 
the early attempts (e.g. Brandenburg & Sokoloff 2002) to 
the most recent (e.g. Courvoisier & Kim 2009; Brandenburg 
et al. 2010; Tobias, Dagon & Marston 2011). Realistically, 
a reliable determination for specific types of stars remains 
a remote possibility. 

With these assumptions on the rotation profile and the 
alpha effect, we solve the ID dynamo equations ^ and 
With Do < then, as expected, the activity wave propa- 
gates equatorwards. The main effect is with the variation 
of TO, larger values of to move the migration nearer to the 
equator. 

With Dq > 0, to = then, again as expected, there is 
polewards migration, centred on mid-latitudes, both with 
j{Q) = 1 and with f{9) nonuniform. When f{9) is nonuni- 
form, the butterfly diagram is concentrated at high lati- 
tudes for modestly supercritical Z?Oj while solutions become 
steady for more supercritical values of Dq. 

With TO = 2 and f{9) = 1, again behaviour is much 
as expected, but quite unexpectedly with f(9) as in Fig. [5] 
and TO = 2,4, the solutions are steady for only marginally 
supercritical values of Dq. These statements are all for solu- 
tions with dipole parity, but steady solutions are also found 
when quadrupole parity is enforced. 



Except when explicitly stated, all the above experi- 
ments were made with a slightly supercritical \Dq\ and 
with a decay term which mimics radial diffusion for the 
value of the ratio /i^^ = h/R = 1/3, where R is the stel- 
lar radius and h the thickness of the convective zone, e.g. 



Baliunas et al. (2006) We deduce that, even with this sim- 
ple dynamo model, behaviour can be remarkably rich and 
broad conclusions cannot be drawn without a quite careful 
exploration of the parameter space and the form of f{9). 

3.2. 2D dynamo models 
3.2.1. Solar rotation law 

Now we turn our attention to 2D models and start by using 
a solar- like rotation law within fractional radii 0.64 < r < 
1.0, as shown in the left hand panel of Fig. [TJ In r > 0.7 
this is an interpolation on MDI data, and is made to match 
uniform rotation at the lower boundary r = = 0.64. 

A dynamo shell with a lower boundary at vq = 0.64 
of the stellar radius is a good assumption for G and K 
type main-sequence stars, but it is not appropriate for a 
subgiant such as the KlIV active component of HR 1099. 
Given the current uncertainty on the evolutionary stage 
of the star, rp = 0.2 is a reasonable assumption (cf., e.g., 
Lanza (2005)). Because extensive numerical simulations to 
cover realistic values for this (and many other) dynamo 
governing parameters are obviously beyond the scope of the 
paper, we restricted our investigation here to extrapolations 
that seem reasonable on the basis of available knowledge. 

Several forms of the functions fi{r),f2{9) are used: 
f2{9) = cos sin™ 0, with to = 0,2,4, and /i(r) takes 
the forms shown in Fig. [31 referred to below as ia — 
0, 1, 2, 3, 6, 7, as denoted in the caption. There are thus 18 
possible forms of a{r, 9). 

A standard value — ilcqR^/vo = 6 x lO'^ was taken. 
The other dynamo number, Cq,, was given a slightly super- 
critical value. The results are summarized in Table 2. 

Only odd parity (dipole-like) solutions were studied. 
Models 1-18 have Ca < 0, the conventional fix to get near- 
surface fields migrating in the solar sense. Models 19-30 
(not all numbers present) have Cq. > 0. Latitude-time dia- 
grams for the pairs (1,19), (6,24), (8,26), and (12,30), which 
have the same values of the parameters ia and to and the 
same |Cq|, are shown in the top four rows of Fig. |4l Here 
and below, "near-surface" refers to a fractional radius of 
ca. 0.93, and "deep" to radius ca. r = rg. 



3.2.2. Meridional circulation 

Here we investigate the effects of an arbitrarily imposed 
meridional circulation, in addition to solar-like differential 
rotation and a-effect, on the butterfly diagram. 

We take a circulation determined by a stream function: 



= Rm-{r — ro)'^{r — 1) sin^ 0cos0, 



(5) 



so that 



9V 



ue 



1 dip 
• sin 9 dr ' 



(6) 
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Fig. 3. The function /i(r) (which determines the radial dependence of the unquenched a-term). Row (a), io 
(b), ia — 2, 3; row (c), ia = 6, 7. 



0.1; 



Table 2. The general features of the solutions with solar-like rotation law. EW, PW, SW denote equatorwards, polewards, 
standing wave respectively for migration of in a near-surface layer. Models 1-18 have Cq < 0, Models 19-30 have 
Ca > 0. Models (1,19), (6,24), (8,26) and (12,30) have identical parameters, except that the sign of Ca has been reversed, 
i.e. the pairs have the same value of \Ca\- 



Model 




m 


Notes 


1 


3 


4 


predominantly EW at low latitudes, weaker PW branch at high latitude 


2 





4 


rather strange butterfly, very weak EW at low lats, no migration at high lats 


3 


2 


4 


low lat EW, no high lat feature 


4 


6 


4 


steady solution 


5 


7 


4 


steady solution 


6 


1 


4 


predominantly EW at low latitudes, very weak PW branch at high latitude 


7 


3 


2 


predominantly EW at low latitudes, weaker PW branch at high latitude 


8 





2 


strong PW branch at high lat, EW at low 


9 


2 


2 


low lat EW, very weak high lat PW 


10 


6 


2 


steady solution 


11 


7 


2 


steady solution 


12 


1 


2 


EW at low lats, weak PW at high lat 


13 


3 





weak EW at low latitudes, much stronger PW branch at high latitude 


14 








strong PW at high latitudes, very weak EW branch at low latitude 


15 


2 





steady solution 


16 


6 





EW at very low lats, no high lat features 


17 


7 





steady solution 


18 


1 





strong PW at high lat, very weak EW at low lat. Solus steady at large Ca 


19 


3 


4 


almost no migration - near SW low lats, with very weak PW drift 


24 


1 


4 


much as Model 19 


26 





2 


EW, extending over most latitudes 


30 


1 


2 


mild EW (near SW), over most latitudes 



where tq = 0.64 corresponds to the base of the dynamo re- 
gion. (Taken literally, in the solar case this implies the cir- 
culation penetrating into an overshoot region, but limited 
experimentation with the circulation restricted to r > — 
0.7 suggests little difference). Here Rm — UoR/ijo, where 
Uq is the maximum value of ug at the surface. This circu- 
lation has a single cell in each hemisphere, with polewards 



flow at the surface if Rm > 0. The streamlines = const 
are shown in Fig. [Sj 

In our models, the a effects and the maximum of 
the angular velocity shear are not spatially separated, 
as in, e.g., flux transpor t models based on the Babcock - 
Leighton paradigm (e.g., Dikpati & Charbonneau (1999)). 
Therefore, the meridional how does not play the crucial 



6 




a 



a 



a 



a 




-90 i 



5.25 

time 



-90 J 

0.201 



0.3 



time 



Fig. 4. Butterfly diagrams for the sub-surface fleld. Left hand cohinin has Ca < 0, right hand column has Ca > 0. 
Rows 1-4 have the quasi-solar rotation law, rows 5 and 6 the quasi-cylindrical law. First row: Models 1 and 19 with 
ia = 3, m = 4. Second row: Models 6 and 24 with = 1, m = 4. Third row: Models 8 and 26 with = 0, m = 2. Fourth 
row: Models 12 and 30 with ia = l,m = 2. Fifth row: Models Ic and 13c, with = 3,m = 2. Bottom row: Models 7c 
and 19c, with = 3, m = 4. 
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1 









90 
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Fig. 5. Streamlines of the circulation defined by With 
Rm > the surface flow is polewards. 



role it has in those kinds of models, but it can still modify 
the migration of the activity waves in a given layer when 
its speed there becomes comparable to the migration speed 
of the wave in the absence of circulation which is mainly 
established by the product \adi}/dr\. 

We first look at Model 26 (with a solar-like rotation). In 
the absence of any meridional circulation, the near-surface 
migration is equatorwards (Tabled, but the deep butter- 
fly diagram has polewards and equatorwards branches, the 
polewards being somewhat stronger, see Fig. [51 top panel. 
In the presence of our one-cell circulation, polewards at 
the surface, the near-surface migration remains equator- 
wards for Rm ^ 40, but the deep equatorwards branch 
is strengthened, see the lower panels of Fig. [Bl For larger 
values of Rm, the near-surface migration acquires a weak 
polewards branch when Rm ^ 100, but the deep migration 
develops two discrete patterns that are close to standing 
waves. Note also that, here and below, we do not claim to 
have explored exhaustively the parameter space, but just 
to have sampled a few, we hope fairly representative, solu- 
tions. Other behaviours may well remain to be found. 



3.2.3. Quasi-cylindrical rotation law 

Here a quasi-cylindrical rotation 



law, 



Covas, Moss & Tavakol (2005) and shown in the right hand 
panel of Fig. 1, is used. We takeCa = -11.7, = l.SxlO'^, 
i.e. significantly supercritical values that might correspond 
approximately to = Sfi©. Exceptionally, Models 13 
and 19 have Cq = -1-11.7 > 0. Results are summarized 
in Table 3 (note that these models are labelled "Ic" etc 
to distinguish them from the models with a solar-type 
rotation law). Butterfly diagrams for the sub-surface 
toroidal field for the pairs of Models (Ic, 13c), and (7c, 
19c) are shown in Fig. U] in the fifth and sixth rows, 
respectively. The models in each of these pairs have the 
same parameters, except that for 13c and 19c, Ca > 0. 

3.2.4. An intermediate rotation law 

Although several stars in our study are quite rapid rotators, 
the cylindrical rotation law discussed in the previous sec- 
tion is not the only possibility for their rotation. It appears 
plausible that a regime intermediate between the solar-like 
and the quasi-cylindrical rotation laws may be more appro- 
priate for such stars. Therefore, we explored the butterfly 
diagram obtained by adopting such a synthetic rotation 
curve, i.e., mixing these laws with weights of about 50% 
(Fig. [7]). Using this il(r, 9), we obtained butterfly diagrams 



S 



■O 
3 




4.251 



4.27 4.28 

time 



Fig. 6. Butterfly diagrams for Model 26 (solar rotation, 
Ca > 0,ia — 0,m = 2) : top - without meridional cir- 
culation (Rm — 0), bottom of the convection zone; mid- 
dle - Rm — 20, sub-surface; bottom - Rm = 20, bottom 
of the convection zone. Solid/broken contours denote posi- 
tive/negative values of respectively. 




Fig. 7. Synthetic rotation curve from combining solar-like 
and quasi-cylindrical rotation laws with 50% weighting. 



for Bcf, in both the deep and shallow parts of the convec- 
tive zone. These are characterized by field migrating from 
low and high latitudes towards some intermediate latitude, 
even without including any meridional circulation (Fig. [5]). 
Neither the solar-like nor the quasi-cylindrical rotation law 
produced such a behaviour. 



3.2.5. Models with a deep convection zone 

Since tq = 0.64 probably gives a CZ that is too shallow for 
HR 1099, we experimented with models having a deeper 
dynamo zone (rg = 0.2), with both solar- like and quasi- 
cylindrical rotation laws and also their equally weighted 
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Fig. 8. Butterfly diagrams for the synthetic rotation law 
(combination of solar-like and quasi-eylindrical) : top panel 
- shallow, lower panel - deep parts of the convective zone. 



combination. The role of a simple meridional circulation, as 
parametrized by Rm, was also considered (cf. Sect. 13.2."^ . 

In general, the deep convective zone appears to favour 
steady solutions when Rm ^ 0, or standing wave (SW) so- 
lutions. The latter may partially be due to the restricted 
latitudinal extent at the bottom of a deeper convective 
zone. This conclus ion agrees with the previous results of 



Moss et al. (2004) 



The solar rotation law sometimes gives low latitude 
equatorwards and high latitude polewards branches simul- 
taneously present near the surface, similar to the models for 
solar-like stars with tq = 0.64 previously discussed. With 
Cq < 0, with the quasi-cylindrical rotation law, deep and 
surface migration was either SW or equatorwards depend- 
ing on the sign of Rm, and becomes steady for large enough 
|i?m|. If Cq. > 0, then solutions with polewards surface mi- 
gration and SW or vacillatory (V) behaviour deep down 
were found for Rm > 0. With Rm < 0, for the same Ca 
surface migration patterns were a combination of SW and 
equatorwards, again with SW in the deep CZ. Again solu- 
tions are steady for large enough |i?TO|. The richest choice 
of possibilities was obtained for the intermediate rotation 
law with Ca and Rm varying (Table |3]). It is notable that 
examples with surface polewards migration were not found. 
These investigations clearly demonstrate that the effects of 
advection can dominate the basic dynamo wave when Rm 
is significantly non-zero. We conclude that none of these 
butterfly diagrams can reproduce the phenomenology ob- 
served in HR 1099 (case I of Sect. 

A more promising result fFig.llOp can be obtained when 
using an a profile which changes sign with depth (see 
Fig. in]) which gives clearly opposed waves near top and 
bottom of the convective zonesij. In this case, we only inves- 



tigated the quasi-cylindrical rotation law for a very limited 
range of dyna mo parameters. This resu lt is consistent with 
the results of Moss & Sokoloff (2007) concerning dynamo 
waves propagating in two separate layers, and allows the 
behaviour sought with appropriate choices of Ca in the two 
layers. The separation o f the layers, that was assumed by 
Moss & Sokoloff (2007) arbitrarily, may be achieved here 
by exploiting the freedom given by the depth of the con- 
vective zone. Of course, this assumes that the deep and 
surface magnetic fields somehow jointly contribute to the 
surface activity manifestations (see Sects. S] and [5]). 



4. Comparison of the 2D dynamo models with the 
observations 

With the solar rotation law, details of migratory patterns 
are sensitive to the spatial dependence of a, and even the 
general rule that the direction of migration is governed by 
the sign of CaC^ appears not always to hold true. The 
most striking case is that of Models 8 and 26 - see the 
third row of Fig. IH Model 8 (left panel) has Cq, < and 
shows a strong polewards branch in addition to an equa- 
torwards low-latitude branch. Reversing the sign of Ca but 
maintaining the same rotation law and spatial structure 
of the a effect, we see only a pronounced equatorwards 
branch, contrary to simple intuitive expectation. The effect 
of switching from Cq < to Cq > is also remarkable and 
not simply intuitive for Models 6 and 24 (on the second row 
of Fig. and for Models 12 and 30 (on the fourth row); 
these differ in the latitudinal localization of the a effect 
(to = 4 and to = 2, respectively). 

With the quasi-cylindrical rotation law, migratory pat- 
terns are quite insensitive to the exact form of a, the only 
significant change occurs when the sign of CaC^j is reversed. 
The rule linking the sign of CqC^j to the migration direc- 
tion appears to hold. Plausibly this is because the spatial 
structure of 17 is much simpler than in the solar-like case, 
and a cannot be 'tweaked' so as to give extra weight to re- 
gions of the envelope with anomalous gradients of i7, unlike 
in the solar case. 

Our general impression from the above analysis can be 
summarized as follows. Dynamo models with rp = 0.64 and 
a fixed sign of a through the CZ can provide polewards mi- 
grating patterns, and observers understand activity in some 
stars as a manifestation of a polewards propagating pattern 
(see cases III and IV in Sect.[2|). Models with a solar-like ro- 
tation law tend to show both equatorwards and polewards 
branches, the latter with an intensity that depends on the 
spatial distribution of the a effect. It is possible to repro- 
duce case IV by, e.g.. Model 8, while the case with a stronger 
equatorwards and a weaker polewards branch, reminissent 
of the behaviour observed in the Sun, can be compared with 
Model 6 (see Fig. m left hand panel in the second row). 
The case with polewards migration only, i.e. case III of our 
classification, is unlikely to be reproduced with a solar-like 
rotation law, but can be reproduced with a cylindrical law 
(cf. e.g. Model 13c with Cq > 0, the right hand panel, fifth 
row of Fig. H]) . This rotation law may be characteristic of 



^ Turbulent convection simulations by Kapyla et al. (2009) do 
indeed show an a effect that changes sign from the bottom to 
the top of the convective domain (cf. the sign of their tensor 



coefficient Ofyy, that corresponds to the scalar coefficient a of our 
simplified mean-field model, in their Figs. 3 and 12). This effect 
also appea rs when a is determined from some tur bulence models 
- see e.g. 
therein. 



Brandenburg & Subramanian (2005) and references 
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Table 3. The general features of the solutions with the quasi-cylindrical rotation law. EW, PW denote equatorwards, 
polewards respectively for migration in a near-surface layer. Note Ca — +11.7 > in Models 13 and 19. 



Model 




m 


Notes 


Ic 


3 


2 


strong EW at low latitudes, no high latitude features 


2c 





2 


strong EW at low latitudes, no high latitude features 


3c 


2 


2 


strong EW at low latitudes, no high latitude features 


4c 


6 


2 


strong EW at low latitudes (slightly more extended than 
above), no high latitude features 


5c 


7 


2 


similar to Model 4c 


6c 


1 


2 


similar to Model 4c 


7c 


3 


4 


similar to Model Ic 


8c 





4 


similar to Model 2c 


9c 


2 


4 


again, low lat EW 


10c 


6 


4 


similar to Model 9c 


He 


7 


4 


as Model 10c 


12c 


1 


4 


as Model 10c 


13c 


3 


2 


Ca > 0. PW at mid-latitudes 


19c 


3 


4 


Ca > 0. PW at low latitudes 



stars rotating significantly faster than the Sun (see, e.g., 
Lanza (2005)| . 

On the other hand, a pattern of spots extending from 
the high to the low latitudes without any clear evidence 
of migration during the activity cycle (case II in Sect. [U, 
may be interpreted by, e.g. Model 24 (see Fig. HI right hand 
panel in the second row) . That model is characterized by a 
stationary pattern spanning a latitude range from the equa- 
tor up to ~ 65° - 70°. Also Model 19 (Fig. H right hand 
panel, first row) shows a similar stationary pattern, but 
more localized in latitude. Model 12 (Fig. 01 left hand panel 
of the fourth row) shows a weak polewards branch that, 
although probably not capable of producing photospheric 
spots, could be detectable through the chromospheric Ca II 
H&K flux modulation. Therefore, it might correspond to 
our case VI because the pattern localized in latitude gives 
rise to an almost constant primary periodicity in the chro- 
mospheric flux modulation, while the migrating branch pro- 
duces a secondary periodicity that varies along the activity 
cycle. 

When interpreting the observations, one should consider 
that the outer contours of the model butterfly diagrams can 
be affected by the turbulent diffusion. Nevertheless, the di- 
rection of migration, which is the relevant information for 
our study, is not modified, as can immediately be seen by 
considering the direction of migration given by the interior 
contours of the diagrams that correspond to increasingly 
stronger fields. In some of the diagrams, such as those of 
Model 6 (on the left second row in Fig. 2]) or Model 12 
(on the left fourth row) we plot only one contour to repre- 
sent the polewards branch, but the reality of that branch is 
confirmed by more detailed investigation. (A single contour 
merely indicates that the field strength in that branch is rel- 
atively small.) Moreover, turbulent diffusion in our model 
is uniform in the outer layers of the convection zone which 
implies that it cannot favour a specific direction of migra- 
tion of the field, i.e. it does not introduce any preference for 
equatorwards or polewards motion of butterfly contours. 

In the case of HR 1099 (case I of Sect. [2]), the point 
is that the image of polewards patterns extracted from 
the observations is quite different from that emerging from 
the theory in all the models with vq — 0.64 introduced in 
Sect. 13.21 irrespective of the adopted rotation law or the in- 
clusion of a meridional circulation. We recall that HR 1099 



displays two patterns that migrate in opposite directions, 
and the migration that begins nearer to the equator is pole- 
wards. 

If there is simultaneous migration polewards and equa- 
torwards through the same latitudes, then it is difficult to 
see how any simple mean-field model can reproduce it. We 
attempt to interpret the two oppositely propagating activ- 
ity patterns as originating from different spatial volumes 
(see Sect. I3.2.5|) . Specifically, in the case of HR 1099, we 
need to consider an a effect that changes sign with radius 
or the presence of two distinct dynamo layers with opposite 
signs of a. The main difference with respect to the models 
with 7'o — 0.64 is that now we consider the contribution 
to surface activity from both the deep and the shallow dy- 
namo layers, while in the other cases we assumed that the 
field pattern at the top of a single dynamo layer directly 
accounts for the observed butterfly diagram. 

We recognize that there are unresolved difficulties with 
this idea. For example, there is the problem of storing 
toroidal fields in relatively shallow superadiabatic lay- 
ers without a too rapid loss through buoyancy insta- 
bilities. This issue has been addressed, inter alia, by 
Brandenburg (2005, 2009). We note in passing that in 
our model the toroidal field strength in the upper region is 
about 20% of that in the lower. This might assist downward 
turbulent pumping to stabilize the field by its reduced mag- 
netic buoyancy force (proportional to B^), possibly acting 
preferentially at the base of a stellar supergranular layer 
analogous to the solar supergranulation (but at a some- 
what greater depth because the convection in a subgiant 
is expected to have larger vertical scales). Turning back 
to the Sun, other difficulties for a dynamo model operat- 
ing in the subsurface shear layer may include an incorrect 
phase difference between poloidal and toroidal fields and 
a too weak toroidal field due to the strong turbulent dif- 
fus ion expected in th ose highly superadiabatic layers, cf. 
Dik pati et al (2002)[ 

Nevertheless, we can attempt to reverse the argument: 
from both general considerations about mean-field models, 
and our particular simulations, we are unable to identify a 
mechanism to produce the simultaneous presence of pole- 
wards and equatorwards migration co-located in latitude, 
other than that discussed above. Thus we conclude that 
either this or a related mechanism does indeed operate, or 
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that the phenomenon is beyond the scope of mean- field the- 
ory. 



Note that [Berdyugina fc Henry (2007) stress the non- 
axisymmetric behaviour of HR 1099 and other active stars, 
whereas we restrict ourselves to axisymmetric models. We 
fully realize that departures from axisymmetry may play 
a role in the phen omena displayed by HR 1099, cf. e.g. 
Moss et al. (2002) [ and seem essenti al to explain the 'flip- 
llop' phenomenon, e.g. Moss (2004) However nonaxisym- 



metric dynamo models contain many additional uncertain- 
ties and the main features of the observations seem to be 
reproduced by a simpler model. 



5. Discussion and conclusions 

Our general conclusion is that it is becoming realistic to 
construct a framework for classification of the very var- 
ied stellar dynamo wave behaviours, and to relate this to 
physical parameters of stellar convective zones. We have 
presented above a very crude and preliminary outline for 
such a template. 

We confirm that the sign oi D = CaCui ^ adVl/dr is 
the main quantity which determines the direction of ac- 
tivity wave propagation, even though considerable finer 
detail can be found in the results, e.g. in the form of two 
branches at low and high latitudes, or of standing wave 
patterns (cf. Sect. [Hand Fig.|H). However, quite unexpect- 
edly, we conclude from the above plots that even if Z? > 
it is quite difficult to excite a pronounced single polewards 
branch in stars with a solar-like rotation law. It follows that 
straightforward considerations based on the sign of dynamo 
number D are inadequate to understand the occurence of 
polewards branches, and a careful examination of 2D mod- 
els is important. In our investigation, we consider only the 
direction of dynamo wave migration and do not attempt 
to fit the observed spot latitudes. The direction of migra- 
tion (or standing wave behaviour) appears to be a robust 
result that does not depend on the details of the butterfly 
wings (which may be influenced by physical processes that 
we have not considered). 

Observations give a hint that one hemisphere of a 
star can contain two oppositely propagating activity waves 
which are pronounced enough to be observable. We find 
that, at least for some stars, rotation curves and spatial dis- 
tribution of the other dynamo governing quantities should 
produce two activity patterns in a hemisphere, of more or 
less comparable intensity. A phase difference between the 
equatorwards and the polewards dynamo waves can exist. 
This idea emerges from looking at the butterfly diagram 
for Model 1 (top left panel of Fig. 2]), where there is a dif- 
ference in the field intensity of the two branches as well 
as a phase shift between the epochs when the low-latitude 
branch reaches the equator and the high-latitude branch 
reaches the highest latitude. 

The fundamental result of our investigation is that the 
spot migration observed in main-sequence stars or the pres- 
ence of standing activity waves can be explained by consid- 
ering the time evolution of the magnetic field at the upper 
boundary of a dynamo shell, even in a very simple dy- 
namo model. The various behaviours can be reproduced by 
changing the spatial structure of the a-effect and/or the 
internal rotation law. Also a meridional fiow may play a 
role. 




0.2 0.4 0.6 0.8 

radius 



1.0 



Fig. 9. Dependence of a on radius at 6* = 7r/4 for the case 
where a changes sign with radius. 
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Fig. 10. Butterfly diagrams resembling that of HR 1099 for 
a deep convective zone with a change of sign for a (Fig. [9]), 
quasi-cylindrical rotation, —Q: sub-surface (upper panel) 
and deep (bottom panel). 



The other important point is that the behaviour of 
HR 1099 cannot be explained by our simple one-layer 
models and the inclusion of a meridional circulation 
does not change this conclusion. On the other hand, a 
tw o-layer model, somewhat similar to that introduced 
by Moss fc Sokoloff (2007)^ can explain the behaviour of 
HR 1099 if the two spots migrating in opposite directions 
are the result of magnetic flux tubes originating in the 
upper and deeper dynamo layers, respectively. However, 
this interpretation requires a significant modification of the 
paradigm used for main-sequence stars. In other words, 
while for main-sequence stars we use a single-layer dynamo 
and assume that the observed spots correspond to the field 
at the upper boundary of the dynamo domain, for the sub- 
giant in HR 1099 (and for subgiant stars in general) we 
must assume the presence of two dynamo layers separated 
by an inactive shell, with both layers contributing to the 
observed spots at the surface. 

From the viewpoint of dynamo theory, we have shown 
that by exploiting the considerable range of freedom in 
choosing ill-known physical quantities, we can produce a 
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Table 4. The variety of butterfly diagrams obtained with the intermediate rotation law and the deep convective zone 
(tq — 0.2). Some intermediate cases, not tabulated, were also found. These behaviours were each found for a variety of 
parameters and, given the fragmented nature of the parameter ranges investigated, we have not here specified parameters, 
and rather just indicate the possible range of behaviours. 



near surface 


deep 


near surface 


deep 


weak PW at mid-latitudes, 
SW at high lats 


SW 


V/SW near poles only 


V/confused 


near SW 


weak EW 


SW, concentrated at high lats 


weak PW, mid-lats 


EW 


almost no fleld present 


EW 


PW 


steady 


steady 


steady 


PW 


EW 


EW 


PW 


EW 


confused butterfly, no migration 


near SW 


EW 


SW 


V/confused, near pole only 


V/confused 


EW and SW 


EW and SW 



wider range of activity wave behaviour than had previously 
been realized - indeed we can find something approach- 
ing the observed range of stellar surface activity waves. Of 
course, theory cannot at present say which of the models 
(if any) correspond even loosely to reality. We just make 
the point that the various observed behaviours are not in- 
compatible with even simple dynamo models. The inherent 
uncertainties of mean field models make it difficult to make 
a stronger or more useful statement than this. It remains to 
be clarified, for example, whether there is a correlation of 
the dynamo wave behaviour with the absolute value of the 
dynamo number, as suggested by our ordering of the cases 
in Sect. [2]by increasing values of Ro. To take just one point 
considering a solar-like rotation, when the dynamo num- 
ber is large (and Ro is small) we generally see two waves 
with comparable magnetic field strength, starting from mid 
latitudes and migrating towards the pole and the equator, 
respectively. On the other hand, when the dynamo num- 
ber is small (and Ro is large) the equatorwards wave has a 
stronger field than the polar wave. Further observational 
evidence and characterization of polewards waves would be 
of special interest and value. Also the simultaneous presence 
of a standing dynamo wave at a fixed latitude together with 
a migrating wave is of interest in interpreting the behaviour 
of some of the Mt. Wilson stars. 

Finally, we have of necessity based this exploratory pa- 
per on a rather small number of relatively well observed 
stars (and of a small subset of all possible dynamo mod- 
els). We thus make the usual plea for more, reliable, data 
in order to substantiate (or disprove) our attempts at sys- 
tematization. 
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